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ABSTRACT : 

The  predicted   dependence   of    sound   transmission    on    the   statistics 
of   the   randomly-rough    interface   between    dissimilar    fluids    has    been 
studied   by   use   of    the  Helmholtz    Integral    (J.    Acoust.    Soc.    Amer .    51, 
IO83-IO9O    (1972)).      The  predictions   have   been   verified  for    radiation 
from  a    helicopter    hovering,    and  slowly   moving,    over    the  sea,    for    fre- 
quencies   to   1000   Hz    for    a   wide   range   of   surface   acoustical   roughnesses, 
R   =   k|tf    [  (C2/C1)  cos   Q-^-cos   ©2]    .      The  roughness   parameters   (7,k,c   and 
©   are   the   rms    height    of   the  surface,    propagation   constant,    speed   of 
propagation   and  angle  with   the   normal,    respectively;    subscript    1   re- 
fers   to  air    and  2   to   the   water. 

When   the   surface    is    mirror -smooth,    the    transmitted   sound  pressure 
is    due   to   the   change    in    divergence   and  change   in    impedance  at    the 
interface.      For    low  roughness,    R  <    1,    the  mean   square   transmitted 
pressure   becomes    decreasingly   coherent   and   is   a    function    of   frequency; 
its   magnitude    is    decreased  by   the   factor    e"R    compared  to   the  perfectly 
smooth   surface.      For    R  ^    1   the    incoherent    component    dominates   and   the 
transmitted  pressure   depends   also  on   the   correlation    length   of    the 
surface   displacements. 

The   transmission   change   of   sound  pressure  as   a   function   of   fre- 
quency   is  presented  for    several    conditions    of  an   SH3-D  helicopter 
hovering  and   flying   over    or    near    an   array   of   microphone   and  sonobuoy 
hydrophones . 
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I.       THEORY 

The  Helmholtz    Integral   has    been   used,    together    with   the  Kirchhoff 
Assumption,    to  predict    the   dependence   of   the   coherent    (phase  cancelled) 
component,    and  the   incoherent   component    of  airborne   sound  transmitted 
through   a   randomly   rough   ocean   surface    .  The  solution    is    given 

simply    in   terms    of  the  rms   height,   O ,    and  the   correlation    length,    L, 

of   the   surface.      The   development    is    very   similar    to   that   used  by   many 

2-7 
others  to  study    the   scatter ing   of  sound  from  a   rough   surface. 

Reference    (1)    describes    various    forms    of   the  problem.      We   sum- 
marize  these   results   by   using   square   brackets    []    to  repeat    the   appro- 
priate  equations    of   that   work.      We  recall    the   cases: 

A)    Infinite  Plane  Waves    Incident    on  a    Smooth  Plane   Interface    in  which 
the  transmitted  pressure   ratio  at    the   surface    is, 

p2s  P2c2/cos92    -  pc   /cosQ 

T    =   ~ —  =    1    +   n    „     ,„„„rs — ^    ~    ~     7,-..^       =    1    +    R  L3J 


pls  P2C2/C°SQ2    +  PxV0080! 


subscripts    1   refer    to  the  medium  of    incidence. 

2   refer    to  the  medium  of  transmission. 

9  =   angle  measured  with   normal 

P  =    density 

c     =  speed  of  sound 

Pls   =  acoustic  pressure,  incident  from  medium  1,  at  the 
surface 

P2s   =  acoustic  pressure,  in  medium  2,  at  the  surface 

R    =  pressure  reflection  ratio  for  plane  waves 

T    =  pressure  transmission  ratio  for  plane  waves 


B)    Aperture-limited  Plane  Waves    Incident    on   a    Smooth   Plane    Interface 
(see   Fig.    1)    for    which   the  Helmholtz   Kirchhoff   solution   for    a    large 
aperature    (AA  »  K)     is    given   by    equation   [6a]: 


p*.*H<pfc-rfe>-B 


L5l 


Y         X 


P2m    AA 


e  dxdy 


[6] 


where 


=  p„    F  sine   K      X  sine   K   Y 
r2m  x  y 


2m 


-i2p1Qk2e 


ik2r2 


4fTr 


AA   cos    0, 


AA    =    4XY 


ik2r2'  "k2r2      -ik    .; 

„,    _    e _   e 2 

if)    =  . =  e 


r    =    i  x    +    i    y   +   i   C,    the   radius    vector    to   the  surface   element 
-L  ^  3 

r  =   distance   from   source   to   origin   of   ensonified  area 

r    '  =   distance   from   source   to  surface   element 

r  =   distance   from   origin   of   ensonified  area    to  receiver 

r     •  =   distance   from   surface  element    to  receiver 


sin   K  X 
SinC   KxX   =         K   x" 


(c  /c   sin  ©   -  sin  ©  )  tan  ©   cos  ©     sin  ©   tan 

=  1  + ; x x •*■   j 


e2  sin2Q3 


c  /c   cos  ©   -  COS  © 


c?/c   COS  ©  -  COS  © 


p     i(k1r1-0)t) 
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,  incident  acoustic  pressure  at  surface  origin 


p_ .  =  acoustic  pressure  at  unit  distance  from  source 
1A 

p  =  acoustic  pressure  at  the  receiver 

dS  =  dxdy/cos  y ,    area  of  a  surface  element 

C  =  C(x,y),  surface  height  at  (x,y) 


See  below 


GEOMETRY  OF  THE 
PROBLEM 


Figure   1.    Geometry    for    Helmholtz -Kir chhof f   theory    for    transmission 
of  a   plane   wave   through   a    randomly    rough   interface   betwee 
two   dissimilar    media. 
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^  a 


n      =    i    (-sin  y   cos   a)+    i    ( -s  in  y    sin  0i)+    i      cos   V,    normal    to 
J-  ^  3 

surface   element 
V      =   angle    the   normal    makes    with    the   z   axis 
OL      =   angle    the  projection    of   the   normal    onto   the   xy-plane   makes 

with   the   x  axis 

-•  2rrf      *  A 

k      =  (i      sin  ©   +    i      cos   ©    ) ,    vector   propagation   constant    of 


incident    field 


•*  2rrf      *  * 

k_    =  (i_    sin  ©_    cos   ©    +    i0    sin   ©0    sin   ©„+    i_    cos    ©„),    vector 

2 


propagation   constant    of   received   field 
©      =   angle   of    incidence,    measured   from  normal    to    interface 
©      =  angle    that    k     makes    with   the   z-axis 
©      =  angle   that   k      makes    with    the   xz-plane 

-♦-._»*  A  A 

K    =   k      -   k      =    i  K      +    i  K      +    i_K 
1  2  lx  2y  3Z 

K      =   k  sin   ©      =   k   sin    ©      cos    © 

K      =   k   sin   ©      sin   © 
y  2  2  3 

K      =   k, cos    ©n    -   k   cos    ©_ 
Z  1  1  ^  z 


2  2     2 

K         =    (K      +   K    ) 
xy         v    x  y ' 

C)    Aperture-limited  Plane  Waves    Incident    on  a    Rough   Plane   Interface 
(Gaussian   PDF   of   heights,    Gaussian   correlation    function)    for    which 
the   Helmholtz-Kirchhoff   solution    is    given   by 

<P2P2*>  -R      2       .      2  ,    v      .      2  m    v    •      2  ttl2      -R  :  R"      ^^  ^    ,,., 

—  =   e        F      sine     K   X  sine      K   Y   +   F     — —  e        2 e  L10J 

2  x  y  AA  n=i   nn! 

(p2m} 

2    2 
where   R    =  K     or      =   acoustical    roughness    for    transmission 
z 

=   k22CT2    (c2/Cl   cos    ©x    -   cos   ©2)2  [8] 

(7    =   rms    height    of   surface 


D)    Radiation   From  a    Point    Source    in  Air    Incident   on   a    Smooth  Plane 
Air -water    Interface,    for    which    the   sound  pressure   transmission   change 
for    the  mirror    surface,    TC^    (referred  to  plg)    is 

\    /  \  "1 

[l4a] 


P2m 


TCM   =    20    log10   —  =    -    20    log^ 
Is 


c    cos    0 

"   77  — ^   )  (2COS0.,    cos0o 

I  c      cosO      J  I 


where        D  =   depth   =  perpendicular    distance   from  receiver    to   interface 

H   =   height   =  perpendicular    distance   from  source   to   interface 

p^      =  pressure  at    Snell   angle  90   depth    D,    for    a   mirror-like   surface, 
r2m  2 

and  for    which   the  change  referenced   to  the   incident   pressure  at   ground 
zero   (p,p)    is   given   by 

tcmz  =  tcm  "  20  logio  cos  ei  Cl4b] 

E)  Radiation  From  a  Point  Source  in  Air,  Incident  on  a  Randomly -rough 
Plane  Air -water  Interface. 

This  last  problem  is,  needless  to  say,  the  most  difficult.   The 
approximate  approach  derived  in  Reference  1  and  used  here  (Section  III- 
C-3)  is  similar  to  that  in  Reference  7.   The  ensonified  area  is  divided 
into  subareas;  then  all  the  coherent  contributions  and,  separately,  all 
the  incoherent  contributions,  are  added.   The  often  incompatible  criteria 
for  each  subarea  are:  1)  the  incident  pressure  is  approximately  constant 
and  the  spatial  variation  of  the  incident  phase  is  negligible  (e.g.  less 
than  X/8)  over  the  sub-area;  2)  the  linear  dimension  of  the  subarea  is 
several  surface  correlation  lengths  in  extent;  3)  the  correlation  length 
is  much  greater  than  the  acoustical  wave  length.   The  desirability  of 
assumption  (1)  in  adapting  the  plane  wave  solution  is  obvious.   Assumpt- 
ion (2)  states  that  the  infinite  integration  that  yields  [lO]  is  accom- 
plished after  "enough"  correlation  distances  so  that  the  spatial  cor- 
relation has  reached  a  "small"  value.   We  will  say  more  about  this  later 
(Section  III-C-3) .   Assumption  (3)  is  the  requirement  for  our  use  of  the 
Kirchhoff  Assumption,  in  which  the  infinite  plane  wave  impedance  relation 

at  the  surface  is  assumed  on  a  point-by-point  basis. 
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For  the  point  source  we  therefore  start  with  Eq.  [lO]  for  the 

transmission  through  an  element  of  area  AA  which  we  now  call  A., 

1 

2     2 
multiply  by  (p   ./p,  .)>  and  generalize  by  writing  the  sum  for  many 

j—*  111  V  I    ~*   J 

such   areas,    (see   Fig.    2). 

-    <P2iP2i*>      r   P2mi      -Ri      2       .       2V      v         .2         „ 

E   = =  S   — e  F.    sine   K    .X.    sine   K    .Y.    + 

2  i      2  l  xi    l  yri 

Plsi  Plsi 


1 


2 

p2    .       F2  "L         -R.    -        R"        -K2    .L2/An 
E      2nl      _^ e      l  ,,        _i^  e      xyi  [15] 

^      A .        —  —  i  nn  ; 

1  Pi  •      x        n_1 
^lsi 

We  will  assign  i  =  1  to  the  terms  along  the  Snell  (ray  acoustic) 
direction  and  for  simplicity  we  drop  the  subscript  1  for  this  most 
important  direction. 

An  approximate  solution  to  the  problem  is  obtained  if  we  assume 
that  there  is  no  cross -coupling  between  coherent  and  incoherent  com- 
ponents from  each  subarea.   Then  our  problem  becomes  one  of  separately 
evaluating  the  terms 


2 

Pic"       PlcJ      P- 


<P2P2*>    <P2s>     Var  {p   } 

— =-= —  =  — =2 +  =2-  [!6] 


Is  ^ls  ^ls 

where   the  two   terms    on   the  rhs    of   Eq.    Ll6]    represent    the  respective 
sums    of   coherent    and   incoherent    terms    on    the  rhs    of    Eq.    [l5J. 

The   first    term   of   Eq.    Ll5l    is   a    sum   of   coherent   components;      the 
relative   coherent    transmission  along    the   Snell    direction    is 

<P2S>2  ^      -R   fm         Dc^os^Y2 
rr~       e 


Pis 


1   +  He. cose'  <2   cos9l  cosS2»    •  tl7] 


The    incoherent    term    is    obtained  by   summing    the   effects    of    the 
sub-areas.      Adding   this    to   the   coherent    term,    and   taking    10    log      , 
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Source 


p> c. 


Snell  Ray  for  A,- 


Receiver 


Geometry  for  point  source  above  a  rough  surface. 

Figure  2.      Geometry    for    Helmholtz-Kir chhoff   theory   applied   to 
essentially  plane    sub-areas    of   a    diverging  wave  as 
it  passes    through   a   randomly   rough   interface  between 
two   dissimilar    media. 
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2     2       2 
tt  A.L   cos  ©0  .  F.  S(R.K   L) 

♦  S  — ^ 2ml   X    2XX*        \  [20] 

X2  r2.  (cose^cosSj.) 

where  the  quantities  F.  and  S(R.,  K  L)  are  functions  of  the  angles  at 
-i  l       v  l    xy  '  v 

the  area  element,  9.,  .  and  ©_.  (see  Fig.  2  and  Reference  1). 

li       2i  ' 

Equation  [20  ]  states  that,  for  very  low  frequencies  (R  «  1),  the 

-R 
summation  (incoherent)  term  will  be  very  small,  e   *■  1,  and  the 

transmission  change  will  be  given  by  the  geometry  and  pc  mismatch, 

(Eq.  Ll4].   This  pressure  change  may  be  positive  or  negative,  i.e. 

a  transmission  gain  or  loss,  respectively.   As  the  frequency  (and 

roughness)  increases  in  the  regime  0  <  R  <  1,  the  coherent  transmission 

term  will  decrease  but  will  still  dominate  the  incoherent.   For  R  —   1, 

the  increasing  incoherent  term  will  begin  to  determine  the  magnitude 

of  TC  in  Eq.  [20] .   This  will  cause  an  increase  in  transmitted  pressure 

until  a  high  enough  frequency  is  reached  (value  depending  on  L/A  )  for 

the  S(R,  K  L)  term  to  decrease  with  frequency- and  to  again  cause  a 

greater  transmission  loss. 


II.   EXPER IMENT 

A  series  of  laboratory  studies  were  undertaken  as  a  prelude  to  an 
experiment  using  a  SH3D  helicopter  as  a  noise  source  over  the  sea.   The 
laboratory  results  and  two  graphs  of  the  ocean  data  have  been  in- 
cluded in  a  paper  titled  "Spectral  Characteristics  of  Sound  Transmission 
Through  the  Rough  Sea  Surface",  by  H.  Medwin,  R.A.  Helbig  and  J.D.  Hagy, 
Jr . ,  which  has  been  submitted  for  publication  in  the  Journal  of  the 
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Acoustical    Society   of  America    (estimated  publication    date   Dec    1972). 
This    technical   report   contains    the  complete  presentation   of  the 
transmission   results    of   the  sea    trials. 

The  point   source   theory   for    a   rough   interface,    Eq.    L20J,    was    tested 
by   experiments   at   sea    using   the  geometry  as    shown   in   Figs.    3  and  4.      The 
"point"   source  was   a    Sikorsky   SH-3D  helicopter,    at   various    heights, 
hovering   or    slowly  moving   on   command  over    or    near,    an  array    of   one 
microphone   and   several   hydrophones,    on   command.      The   signal    from   the 
receivers   was    carried  by   cable   to  the   nearby   Floating    Instrument 
Platform,    FLIP,    operated  by   the  staff   of   the   Scripps    Institution    of 
Oceanography.      SCRIPPS   scientists    made  measurements    of    the   sea    surface 
height   spectrum,    concurrently  with   this   experiment. 

Two  series    of   experiments   were   conducted.       FLIPEX   I    on   28,    29  May 
1970  and  FLIPEX  II   on    17,    l8  August    1970  both  at   a    deep   sea    location, 
45  miles   NW  of   San   Diego.      Only    the   results    of   FLIPEX   II   have  been 
processed. 

A.      Experimental   Details 

The   SH3-D  helicopter,    viewed  as   a    sound   source,    is    characterized 
by   a    large  number    of  strong   harmonics    of   the   very    low   frequency   blade 
passage   tone   of   the  rotor,    harmonics    of   the  high   speed  counter -rotation 
tail   propellor,    and  a   broad  band  aerodynamic   and  machinery    noise   spec- 
trum,   in  which   the  harmonics   are  partially    immersed.      Only   helicopter 
heights    of  300    ft   or    greater    are  reported  here;    hovers   at    lower    levels 
disturbed  the   sea   surface  and  precluded   comparison  with   theory. 

The   hydrophones   were  parts    of  modified  AN/SSQ-57   sonobuoys.      The 
two  modifications    consisted   of:    a)    disabling   the   RF  modulator,    which 
is    normally   used   to  radio   the   signal    from   the   buoy;    and  b)    changing   the 
length   of   the  hydrophone   support   so   that    it   hung   at    depth   0    (or    10   ft), 
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\                              AIR 
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D^20 

_.__\    RECEIVER 

0, 


Figure  3.   Geometry  of  transmission  into  water  from  point 
source  in  air.   The  Snell  ray  is  indicated. 
(Not  to  scale. ) 
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SIDE  VIEW  OF  EXPERIMENT 
ARRANGEMENT. 


Figure   4.    Experimental  arrangement   for    ocean   experiment        Horn 

source   hovered  over,    or    moved   slowly  past    either    Cluster 
or    Cluster    II.      For    FLIPEX   II,    reported  here,    the   clusters 
were   separated  by    150   ft.      Ocean   and  acoustic    data    were 
recorded  at   Scripps    Institutions   Research  Platform   FLIP, 
300   ft   upwind   from  Cluster    I. 
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20,    40,    95,    200    or    300   ft    below   the   surface.      The   final   configuration 
of  microphone   was    simply   a   hydrophone  placed  on   the  bottom   of  a   4   ft 
diameter   plastic   Kiddy   Pool.      The  colorful   Kiddy   Pool   served  not    only 
as    the    floating   holder    for    the  microphone  but   also  as   an   easily    identi- 
fiable  reference  position   by   which   to   guide   the   helicopter .      The  modi- 
fied  sonobuoy   consisted  of   hydrophone,    pre-amplif ier ,    audio  amplifier 
and  transformer    powered  by   a    salt-water    battery.      Three   sonobuoys   were 
held  together    in   a    cluster.      There   were   two  clusters   as   shown    in   Fig.    4, 
one   300   ft    from   FLIP,    the    other    either    adjacent    to   or    separated  from 
the   first   by  an   additional    in    line,    150   ft.      The  microphone  and  hydro- 
phones   were  hard-wired  to  a    tape  recorder    in   the    instrument    room   on 
FLIP.      The   electrical   cable  was    floated   on   one   inch  polypropylene  line. 
Additional    details    of  the   mechanical,    electrical   and  acoustical   design 
and   tests   will   be   found   in   Ref.    11. 

The  signal    from   the   sonobuoys   was    FM  recorded   on  an   8   channel 
PI   6208   magnetic    tape  recorder,    on   board  FLIP.      One   channel   was   re- 
served  for    voice   annotation. 

The   data   recorded  during   FLIPEX   II   consisted  of   calibrated,    time- 
varying,    voltages    from  hydrophones   at    5   depths   and  a   surface   microphone 
during   helicopter    hovers   at   heights    40,    600  and  750   ft   and   5   knot   flybys 
and   fly-overs   at    150,    300,    600  and  750    ft. 

Under    the   direction   of  Dr.    Russ    Davis    of  Scr ipps ,    time   series   wave 
height    variations    were  recorded  during   both  FLIPEX   I   and   FLIPEX   II. 
Whenever   possible,    such   data   were  taken   simultaneously  with  helicopter 
noise    data.      During   FLIPEX   II,    four    wave  probes   were  used.      Surface 
wave    frequency    spectra   were   generated  by    transforming   the  wave  height 
variation   time   records    for    each   device.      Then   mean   square  wave  heights 
for    several   wave  period   ranges   were    obtained  by    integration    of   the 
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frequency  spectra  over  appropriate  frequency  ranges.   This  processing 
was  carried  out  at  Scripps. 

B.   Data  Analysis 

1.   Analogue  Analysis 

The  frequency  range  100  Hz  to  1000  Hz  corresponded  to  an 
acoustical  roughness  range,  R,  of  approximately  0.05  to  4.0  for  nor- 
mally incident  sound  and  for  the  typica  1  sea  surface  encountered 
during  the  experiment.   Sixteen  frequencies  in  this  frequency  range 
were  selected  for  investigating  the  effect  of  a  broad  range  of  surface 
roughness  on  transmission  and  for  comparison  with  theory. 

The  analysis  utilized  a  Hewlett  Packard  Model  3590A  wave 
analyzer  with  a  General  Radio  Type  1163A  frequency  synthesizer  used  as 
an  external  local  oscillator  for  the  wave  analyzer.   The  logarithmic 
output  of  the  wave  analyzer  was  used  as  the  vertical  input  to  a  Var ian 
Model  F-100  X-Y  recorder.   This  arrangement  was  particularly  useful  for 
analyzing  and  displaying  the  helicopter  passes  since  the  resulting  plot 
showed  the  ambient  noise  levels  before  and  after  each  pass  as  well  as 
the  total  noise  during  each  pass.   A  series  of  plots  made  for  several 
frequencies  using  a  typical  helicopter  pass  appears  in  Fig.  5.   The 
large  time  variation  in  rms  level  even  with  the  slow  response  time  of 
the  wave  analyzer  is  evident.   Signal  voltage  levels  for  the  helicopter 
passes  were  determined  by  taking  an  eyeball  average  of  the  time  vari- 
ations across  the  peak  of  the  pass. 

Ambient  noise  voltage  levels  between  passes  were  determined  in  the 
same  way.   Since  the  plots  represented  rms  levels  in  dB,  this  eyeball 
average  of  dB  levels  was  something  less  than  the  true  average  rms  voltage 
level.   The  error  here  was  on  the  order  of  0.5  dB.   Since  the  type  of 
time  variation  was  about  the  same  for  the  surface  microphone  as  it  was 
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Figure   5.      Strip   chart   records    of   sixty    second  samples    of 
noise   received    in    six    10  Hz   bands    as   helicopter 
passed  over    the  ^0   ft    deep   hydrophone. 
18 


for  the  submerged  hydrophones,  this  error  canceled  itself  out  when 
taking  the  difference  to  compute  transmission  change.   To  obtain  heli- 
copter signal  voltage  levels,  the  ambient  noise  voltage  levels  were 
subtracted  off  by  using  the  white  noise  assumption.   Signal  sound 
pressure  levels  at  the  hydrophone  were  then  determined  by  using  the 
calibrated  sensitivities  of  the  transducers. 

The  Kirchhoff  assumption  was  made  in  order  to  convert  from  ob- 
served microphone  SPL  to  incident  sound  SPL .   That  is,  it  was  assumed 
that  pressure  doubling  had  taken  place  so  that  the  incident  SPL  was 
obtained  by  subtracting  6  dB  from  the  surface  SPL. 
2 .   Digital  Analysis 

The  major  data  reduction  program  for  the  transmission  study 
was  performed  digitally. 

Figure  6  is  a  general  self-explanatory  flow  chart  of  the 
analysis.   In  addition  to  the  flow  chart,  the  following  statements  will 
help  to  characterize  the  analysis: 

The  sampling  frequency  was  2000/sec.   Four  channels,  three 
hydrophones  and  one  microphone,  were  sampled  simultaneously.   For  each 
flyby  or  flyover  experiment,  data  were  digitized  for  peak  sound  levels 
when  the  helicopter  flew  directly  over  the  mic  and/ or  hydrophones. 
Additional  data  were  taken  in  order  to  study  changing  sound  levels 
within  a  time  span  which  began  approximately  50  sec  before  the  helo 
approached  the  overhead  position  and  ended  approximately  50  sees  after 
the  overhead  position. 

Ambient  noise  levels,  digitized  for  corresponding  hydrophone 
depths,  were  taken  at  least  one  half  hour  preceding  or  succeeding  each 
experiment.   Most  ambient  signals  were  taken  within  minutes  of  an 
experiment  and  spanned  a  time  period  of  approximately  15  seconds. 
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After  proper  formatting  of  the  tape,  processing  of  the  data 
was  begun  using  a  system  of  time  seiies  programs    originally  developed 
by  the  University  of  British  Columbia.   The  first  program  in  the  series, 
program  FTOR ,  employs  the  FAST  FOURIER  TRANSFORM  to  compute  Fourier 
coefficients  from  selected  channels  of  the  digital  tape. 

Program  FTOR  was  run  using  varying  block  lengths  of  data  in 
order  to  arrive  at  an  adequate  frequency  resolution  or  bandwidth.   After 
some  experimentation  using  block  lengths  ranging  from  2  sees  to  1/8  sec, 
the  decision  was  made  to  use  blocks  of  1024  samples  (a  duration  of  ap- 
proximately 1/2  sec),  giving  a  frequency  resolution  of  1.95  Hz. 

Contiguous  blocks  of  2    data  points  were  transformed  for  each 
experiment  and  written  on  a  "Fourier  Coefficient  Tape".   This  tape  was 
then  read  by  program  FCPLOT  of  the  time  series  programs.   The  FCPLOT 
program  provides  the  capability  of  averaging  over  only  selected  blocks, 
e.g.,  over  blocks  spanning  10  second  periods,  or  over  all  blocks  as 
was  the  case  for  ambient  data  .   This  program  was  revised  to  write  power 
spectral  densities  for  given  time  spans  on  magnetic  tape  for  further 
analyses . 

From  this  tape  of  power  spectral  densities  created  by  program 
FCPLOT,  blocks  of  data  corresponding  in  time  were  further  averaged  for 
different  runs  of  a  given  experiment.   For  example,  peak  data  were 
averaged  for  several  runs  and  several  ambients  were  averaged.   Data 
from  different  runs  were  weighted  according  to  the  time  span  they  re- 
presented.  The  averaged  power  spectral  densities,  along  with  their 
associated  ambient  levels,  were  read  by  program  TC,  a  program  which 
calculates  sound  pressure  levels  and  sound  transmission  change  and  plots 
them  against  frequency  and  roughness,  respectively. 

Ambient  noise  was  subtracted  from  helo-plus -ambient  levels. 

Points  were  labelled  "bad"  if  ( helo-plus -ambient ) -(ambient )^  0,  or  if  any 
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SPL  differed  by  20  dB  or  more  compared  to  the  SPL  at  a  contiguous 
frequency. 

An  average  TC  was  calculated  for  every  9  contiguous  frequencies 
(band  of  9.75  Hz).  This  average  was  considered  "good"  if  5  points  out 
of  the  9  were  considered  "good". 

Because  of  many  "bad"  points  caused  by  high  ambient  levels  for 
f  <  100  Hz,  the  "low  frequency  SPL"  was  calculated  for  the  frequency 
band  100  to  200  Hz;  the  "overall  SPL"  was  calculated  for  the  frequency 
band  100  to  1000  Hz. 

C.   Comparisons  of  Experiments  with  Theory 

During  the  course  of  the  acoustical  study,  ocean  height  data  were 
taken  simultaneously  at  four  probes,  at  a  sample  rate  of  eight  times 
per  second,   The  motions  of  FLIP  were  compensated  for.   The  data  for 
each  probe  were  divided  into  four  minute  blocks,  the  spectra  were 
computed,  and  then  the  spectra  were  averaged  across  these  blocks.   The 
wave  height  uncertainty  for  each  of  the  runs  is  probably  less  than  5%. 
The  wave  height  frequency  spectrum  for  each  probe  was  integrated  to 
yield  the  mean  squared  height.   The  values  of  the  mean  squared  height 
at  the  four  probes  were  averaged  and  this  average  was  used  for  the 
theoretical  calculation  of  sound  transmission  changes.   A  calculation 
such  as  described  was  carried  out  for  the  ocean  wave  system  twice  on 
17  August  and  twice  on  18  August  to  yield  values  of  rms  heights,  0 : 
15.5  cm,  14.9  cm,  13.6  cm  and  14.5  cm. 

The  general  comparison  with  theoretical  predictions  is  comprised 
of  several  single  parametrical  checks: 

1 .   Variation  of  Transmission  with  Depth/Height  Ratio  for  Normal 
Incidence  at  a  Very  Smooth  Surface  R  «  1 

When  the  surface  is  relatively  smooth,  R  «  1,  Eq.  [20J  tells 

us  that  only  the  coherent  component  of  the  transmission  change  is 
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significant.   Further,  if  the  surface  is  very  smooth,  R  «  1,  TC 
approaches  TC   (Eq.  l4a)  and  the  transmission  change  will  be  due  only 
to  the  geometry  and  the  pc  mismatch.   The  test  is  accomplished  for  the 
ocean  experiment  by  calculating  the  average  experimental  TC  for  the 
low  frequency  band  100-200  Hz.   For  this  band  the  surface  acoustical 
roughness  for  normal  incidence  sound  (fly-over  or  hover)  was  R  =  0.15 

at  200  Hz  and  one-quarter  that  value  at  100  Hz.   The  decrease  in  co- 

-R 
herent  transmitted  signal  due  to  surface  roughness,  which  goes  as  e   , 

would  therefore  be  expected  to  be  at  most  0.68  dB  at  the  high  fre- 
quency end  of  the  band  and  0 . 16  dB  at  the  low  frequency  end  (100  Hz). 
This  average  expected  loss  of  approximately  0.4  dB  would  not  be  dis- 
tinguishable in  this  experiment.   On  the  basis  of  the  above  argument, 
all  observations  of  transmission  change  for  the  low  frequency  band, 
100-200  Hz,  are  identified  as  occurring  for  R  «  1  and  are  assumed  to 
be  equivalent  to  the  values  that  would  have  been  obtained  for  a  mirror - 
like  surface. 

Figure  7  shows  the  transmission  change  at  normal  incidence  pre- 
dicted by  mirror  surface  theory  (R  «  1)  compared  with  experimental  data 
taken  in  hovers  and  flyovers,  at  a  height  of  600  ft,  for  a  decade  of 
ratios  of  hydrophone  depth,  D,  to  source  height,  H.   Only  the  low  fre- 
quency band  of  sound  frequencies  have  been  used  for  this  test.   A 
positive  value  of  TC  represents  an  SPL  at  the  hydrophone  that  is  greater 
than  the  incident  SPL  at  the  surface.   The  average  disagreement  is 
approximately  2  dB.   This  small  discrepancy  could  have  been  due  to  one 
or  more  of  the  following:  the  helicopter  was  not  a  proper  point  source: 
there  was  a  difference  in  the  absolute  calibration  of  microphone  and 
hydrophones;  the  pressure  was  somewhat  less  than  doubled  at  the  slightly 
rough  surface. 
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2 .   Variation  of  Transmission  with  Angle  of  Incidence  at  a  Very 
Smooth  Surface  (R  «  1)  for  a  Constant  Depth/Height  Ratio 
Figure  8  presents  low  frequency  data  taken  during  fly-overs. 
Each  datum  point  represents  average  over  eight  or  16  blocks  (each 
block  is  %  second  long)  for  cases  in  which  the  aircraft  was  either 
approaching  (at  5  knots)  or  receeding.   The  geometry  and  the  defi- 
nition of  9  n  are  shown  in  Fig.  3«   The  assumption  has  been  made 
(as  above)  that  low  frequency  data  for  which  (R  «  1)  should  be  com- 
pared with  theoretical  curves  for  the  mirror  surface  TC  .   The  com- 

M 

parison    is    again    favorable   for    the   two   D/H  ratios,    O.O67   and  0.5,    for 
which    data   were   available.      The   theoretical   curve   for    D/H   =10,    is 
drawn    for    interest. 

3-      Variation   of  Transmission   with   Sound   Frequency   and  Surface 

Acoustical   Roughness;    Source    over    Receiver 

As    the   sound   frequency    increases,    Eq.    l20J    predicts    that    the 
coherent    term  will    dominate  briefly   and   then   for   R   >    1    the    incoherent 
term  will    control    the   total    transmission    change.       In   order    to   calculate 
the   expected   transmission,    Eq.    [20]    requires    knowledge    of   the   surface 
height    correlation    length  L.      This    can   be   estimated  from   the   known   wind 
speed,    by   using   the   results    of   Cox   and  Munk   to   calculate   first    the   rms 
slope      and   then   using    the  assumption   of  a   Gaussian   Correlation   Length, 
L,    to  calculate    its    value;    this    simple  calculation   yields    the  approxi- 
mate  correlation    length,    L  ~   150   cm   for    the   ocean    during    FLIPEX    II. 

A  second   decision   must   be   made   before   Eq.    C20]    can   be   used   in    the 
case   of  a   point   source;    one   must    settle   on    the  appropriate   sub-areas 
AA.  •      In   principle   they    should  be   chosen   so   that    the    incident    sound 
wave   has    uniform  phase   and  amplitude   over    the   subareas,    and  the   cor- 
relation  has    dropped   to   effectively   zero  at    the  boundary    of   the   sub- 
area.      Figure   9   shows    the   result    of   a    computer    study   of    the   effect 
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of  a)  changing  the  correlation  lengths,  L,  and  b)  changing  the  width 
of  the  circular  rings  AL  which  define  the  sub-area.   The  calculation 
is  for  a  source  at  600  ft  and  a  receiver  at  20  ft  for  a  sea  of  rms 
height,  (J    =    13.  cm  and  correlation  length,  L  ™  150.  cm  (=  5  ft).   It 
is  observed  that  only  the  higher  frequencies  (R  >  1)  are  controlled  by 
these  changes  (that  is,  only  the  incoherent  component  is  modified). 
Changing  the  correlation  length  f r om  L  =  2.5  ft  to  5  ft  to  10  ft  with 
AL  =  3L  has  the  major  effect  on  the  transmission;  changing  the  size  of 
the  sub-area  rings  from  AL  =  L  to  AL  =  3L  or  AL  =  5L ,  for  the  appro- 
priate L  =  5  ft  has  a  smaller  effect  on  the  predicted  transmission  change, 
TC.   Theoretical  calculations  of  TC  to  be  presented,  below,  were  based 
on  L  =  5  ft  and  AL  =  5L  which  fulfills  assumption  (2)  of  Section  I-E. 

It  is  not  possible  to  simultaneously  fulfill  assumptions  (1) 
and  (3)  of  Section  I-E  for  all  the  wave  lengths  studied  in  this  exper- 
iment.  Figure  10  presents  the  experimental  data,  analyzed  by  digital 
computer  (solid  line),  and  analogue  playback,  (circles).   The  two 
methods  of  analysis  are  in  tolerable  agreement.   The  low  frequency  ex- 
perimental values  of  TC  greater  than  +6  dB  are  clearly  due  to  ambient 
noise  underwater.   The  theoretical  predictions  (dashed  line)  appear  to 
have  approximately  the  correct  dependence  on  roughness  (or  frequency 
squared).   However  theory  appears  to  underestimate  the  amount  of  sound 
transmitted  by  roughly  5  dB  for  both  the  hydrophone  at  20  ft  depth  and 
the  one  at  300  ft.   Perhaps  2  to  3  dB  of  this  discrepancy  is  the  same 
as  observed  for  R  «  1  (Figure  7  )  and  the  explanation  of  section  (1) 
could  be  repeated.   However,  it  is  probable  also  that  the  simple  theory 
for  adding  sub-areas  that  was  developed  in  reference  1  underestimates 
the  contributions  of  the  incoherent  sound  in  the  region  R  >  1 . 
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The  transmission  change  for  normal  incidence  sound  is  seen  to 
be  only  slightly  dependent  on  acoustical  roughness  (frequency  squared). 
For  example  the  experimental  curves  could  be  crudely  summarized  as 
showing  a  constant  TC  +  2  dB.   It  is  this  small  spectral  variation  for 
near -overhead  sources  that  has  permitted  some  researchers  to  ignore 
the  frequency  dependence,  completely,  and  to  assume  that  the  mirror 
value  TC^  Eq.  [l4l  as  valid  for  all  frequencies.   This  fortuitous 
situation  is  due,  in  fact,  to  large  incoherent  contributions  from  non- 
Snell  directions  which  compensate  for  the  decrease  in  transmission 
of  the  coherent  component,  as  the  frequency  is  increased. 

The  jaggedness  in  the  experimental  data  is,  of  course,  due  to 
inadequate  averaging  time  being  available.   Much  of  the  data  represent 
averaging  over  10  to  20  seconds,  whereas  10  to  20  minutes  would  be 
necessary  to  reduce  adequately  the  statistical  variations. 

Additional  experimental  curves  of  transmission  change  for  the 
helicopter  over  the  hydrophone  are  given  in  Figures  11-25.   Figures 
11-13  represent  hovers  nominally  over  the  hydrophones  at  20,  40  and 
300  ft  respectively.   In  fact,  the  helicopter  pilot  had  difficulty 
seeing  the  float  marker  and  maintaining  position  during  the  hover  at 
600  ft  altitude  and  unintentional  off-sets  of  as  much  as  100  ft,  may 
have  been  present.   Our  study  of  the  playback  gave  evidence  of  the 
helicopter  drift  and  permitted  the  extraction  of  the  24  second  analysis 

presented  in  Figures  11-13  which  we  believe  represents  the  true  over - 

-R 
head  position.   The  coherent,  e   ,  dependence  for  0.2  <  R  0.5  is  fairly 

clear  in  these  cases  because  of  longer  averaging  time;  the  straight 

-R 
line  starting  at  TC  =  0  shows  the  predicted  slope  of  an  e    dependence. 

The  data  for  R  <  0.2  are  unreliable  because  of  strong  and  variable 

ambient  noise  levels.   The  hover  analysis  yields  virtually  the  same 

TC  data  as  the  fly-over  analysis  shown  in  Figure  10. 
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Figures  14  and  15  represent  5  sec  and  3  sec  analyses  for  source 
at  600  ft  over  the  hydrophone  at  20  ft.   The  150  ft  offset  of  the 
source  relative  to  the  microphone  would  have  a  small  effect  on  the 
incident  levels;  an  inverse-square  correction  has  been  made  for  this. 
The  two  curves  show  the  type  of  variation  of  TC  that  can  take  place 
due  to  a  different  sea  surface  at  the  time  of  analysis  and  due  to  an 
inadequate  time  of  analysis.   Some  of  the  jaggedness  is  eliminated 
when  the  8  sec  average  is  taken;  this  average  was  presented  in  Figure  10 
Figure  16  represents  an  8  second  average  for  the  overhead  source  at 
600  ft  and  receiver  at  40  ft;  it  should  be  compared  with  the  curves 
on  Figure  10  from  which  it  was  omitted  to  avoid  crowding. 

Figure  17  shows  the  transmission  change  for  a  pass  directly 
over  Cluster  II  containing  microphone  and  95  ft  hydrophone.   It  too 
can  be  compared  with  curves  on  Figures  10,  from  which  it  was  excluded 
to  avoid  crowding. 

4.   Variation  of  Transmission  with  Sound  Frequency  and  Surface 

Acoustical  Roughness;  Source  Offset  From  Receiver 

when  the  source  is  not  directly  over  the  receiving  hydrophone, 
our  theory  again  gives  a  fairly  good  prediction  of  experiment.   Figure 
18  shows  the  comparison  of  theory  and  experiment  for  an  offset  of 
approximately  2  50  ft,  resulting  in  an  angle  from  hydrophone,  ©_0  of 
85   for  the  20  ft  receiver  and  40   for  the  300  ft  deep  receiver.   In 
the  case  of  the  20  ft  hydrophone,  simple  mirror  theory  would  predict  a 
transmission  loss  of  more  than  15  dB  for  all  frequencies.   The  theory 
of  Eq.  C20]  comes  closer  to  the  mark  because  it  accounts  for  the  in- 
coherent contributions  which  become  important  at  R  >  1. 

The  region  R  <  0.5  is  of  particular  interest  because  of  the  so- 
called  "evanescent",  "inhomogeneous",  or  "lateral"  wave  that 

13-14 
theorists      predict  would  be  refracted  through  a  smooth  interface. 
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Figure  11.   Transmission  change  in  SPL  as  a  function  of  surface  acoustical  roughness 
for  helicopter  hover  at  600  ft  over  hydrophone  at  20  ft  depth. 
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Figure  12.   Same  as  Figure  11  but  with  hydrophone  at  40  ft  depth. 
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Figure    13.       Same  as    Figure   11   but   with   hydrophone   at    300    ft    depth. 
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Figure    1^+.    Transmission    change   as   a    function    of   surface   acoustical    roughness    for    a 
helicopter    fly-by   at    600    ft    over    hydrophone   at    20    ft    depth. 
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Figure  15.  Same  as  Figure  14  but  a  different  fly-by. 
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Figure   16.    Average   transmission    change   as    a    function    of   surface  acoustical    roughness 
for    two  helicopter    fly-bys   at   600    ft    over    hydrophone   at   40   ft    depth. 
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Figure    17.    Transmission    change   as    a    function    of    surface   acoustical    roughness    for    a 
helicopter    fly-by   at    600   ft    over    hydrophone   at    95   ft    depth. 
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Our  theory  has  not  considered  this  component.   The  fact  that  the  low 
frequency  sound  levels  underwater  are  higher  than  predicted  by  our 
theory  (see  also  Fig.  8  with  D/H  0.067)  may  be  evidence  that  a 
later ial  wave  exists  for  rough  surfaces,  too;  unfortunately,  it  may 
also  be  evidence  of  the  strong  low  frequency  ambient  noise  that  was 
present  during  the  experiment  and/or  refraction  to  the  300  ft  hydro- 
phone due  to  a  strong  (9C  )  negative  temperature  gradient  that  existed 
between  60  ft  and  2Zf0  ft. 

Figures  19-21  show  transmission  changes  as  a  function  of 
roughness  for  a  helicopter  pass  directly  over  the  microphone  but 
offset  by  150  ft  from  the  string  of  three  hydrophones  of  Cluster  I  at 
depths  20,  40 ,  300  ft.   The  helicopter  was  moving  perpendicular  to 
the  line  of  Clusters  I  and  II,  when  these  offset  data  were  collected. 

The  next  series  of  graphs,  Figures  22-25,  show  spectra  of  TC 
at  various  times  during  5  knot  fly -overs  at  600  ft  altitude,  using 
hydrophones  at  depths  20,  40  and  300  ft.   Now  the  helicopter  moves 
perpendicular  to  the  line  of  Clusters  I  and  II  so  as  to  pass  over 
the  hydrophones.   In  Figure  22  we  see  three  curves  for  the  20  ft 
hydrophone  representing  (from  bottom  to  top):  source  50  sec  (410  ft) 
offset  and  approaching;  source  30  sees  (250  ft)  offset  and  approaching; 
source  overhead. 

Figures  23  shows  the  spectra  of  TC  for  the  receding  helicopter 
as  it  continues  to  fly  slowly  along  the  path  perpendicular  to  the  line 
of  the  microphone -hydrophone  clusters:  passing  over  the  20  ft  hydrophone 
in  the  top  curve;  40  sec  (33°  ft)  past  the  hydrophone  cluster  in  the 
middle  curve;  and  50  sec  (410  ft)  past  in  the  bottom  curve. 
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Figure  19.  Transmission  change  as  a  function  of  surface  acoustical  roughness  for 
helicopter  offset  by  150  ft  from  20  ft  deep  hydrophone.  See  text  for 
details.      m 
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Figure  20.  Same  as  Figure  19  but  with  hydrophone  at  40  ft  depth.   See  text  for  details 
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Figure  21.  Same  as  Figure  19  but  with  hydrophone  at  300  ft  depth.   See  text  for  details 
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Figure  22.  Transmission  change  as  a  function  of  surface  acoustical 

roughness  for  H  =  600  ft,  D  =  20  ft  at  different  positions 
of  fly by.   See  text  for  details. 
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Figure  23.   Same  as  Figure  22;  see  text  for  details. 

45 


o_ 

CM 

O  — 

iHa/-^ 

'o- 

vrA 

-~x     , 

y 

Q 

D       ^ 

A/^l 

0 

2. 

0 
ROUGHN 

ESS 

4. 

0 

/ 

/ 
/ 

O- 

jT-s^ 

/ 

1 

^                  ^*- 
^ 

\         /          '^^ 

i 

o. 

CM 

Figure   24.    Same  as    Figure   22,    but   hydrophone 

is    at    40    ft    depth;    see 

text    for    details. 

46 

■ 

Figure  25.    Same  as  Figure  22  but  hydrophone  is  at  300  ft  depth;  see 
text  for  details.      #- 


Figure   24    shows   TC  under    the   same   conditions   as    described   in 
the  previous    figure   except    that   hydrophone    is    40   ft    below   surface. 
The   top   curve   is    for    the   overhead  helicopter;    the   middle   curve   for 
the   source  receding      and  at    165   ft    from  overhead;    the  bottom 

curve    for    the   source     receding    and  at    250    ft    from   overhead. 

Figure   25  presents   TC  for    the   same   conditions    of   flyby   except 
that    the  hydrophone   is    300   ft   below   the  surface.      The   top   curve   is 
for    the   helicopter    overhead  position;    the  middle   curve   is    for    the 
receding      helicopter    now  at   a    distance    of   250    ft    from   overhead;    the 
bottom   curve    is    for    the     receding    helicopter    at   4l0   ft    from   overhead. 
5.      Other    Parameters 

Temperature   refraction   was    an   unavoidable   variable   during    the 
experiments    reported  here.      Figure   26   shows    the   BT  plot    typical    during 
FLIPEX   II.      It    is    believed  that    the  sound  pressures   at    the  20   ft 
hydrophone   were   unaffected  by   the   existing   gradient.      Nor    would  the 
results   at   normal    incidence  be  affected  for    the  hydrophones   at   any 
depth.      On    the   other    hand   the  TC  curves   at    non-normal    incidence   for 
hydrophones    below   20    ft   are   only   approximate   because  we   have  not 
corrected   the   data    for    the   temperature  refraction   effect. 

III.       CONCLUSIONS 
The   Helmholtz-Kirchhoff   theory,    applied  to   transmission    from   air 
into  water,    provides    useful   predictions    of   the   dependence   of   trans- 
mitted sound  pressure   level   on    frequency,    surface   roughness   and  source 
and  receiver   positions.      For    low   roughness,    R  <    1,    transmission    is 
essentially   the   same   as    for    a    smooth   surface,    but    decreased  by    the 
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-R 
factor    e      .      For    larger    acoustical   roughnesses,    R  >    1,    it    is   possible 

to  make  predictions    that   come   close   to  reality  provided   the   ensonified 

surface    is    judiciously    divided   into  sub-areas.      The   experimental    data 

are  not   able   to  clearly      demonstrate   the  presence   of  an    important 

lateral   wave   component    for    low   roughness    surfaces. 

Greater    accuracy    for    the  rough   surface   transmission  problem  will 

require   a    theoretical   approach    that    takes    into  account    the   lateral 

wave   in    the  case   of  a    rough   surface  and  which   considers    interactions 

between    sub-areas    or    deals    with    the  affects    of  diverging  waves    from 

the  beginning  of   the   development. 
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